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We consider the Einstein equation with first order (semiclassical) quantum 

' corrections. Although the quantum corrections contain up to fourth order deriva- 
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equations which contain derivatives no higher than second order. We obtain the 
reduced equations for a range of stress-energy tensors. These reduced equations 



are suitable for numerical solution, are expected to contain fewer numerical in- 
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stabilities than the original fourth order equations, and yield only physically rele- 
vant solutions. We give analytic and numerical solutions or reduced equations for 
particular examples, including Friedmann-Lemaitre universes with cosmological 
constant, a spherical body of constant density, and more general conformally flat 
metrics. 
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1. Introduction 

Quantum corrections to general relativity are expected to be important in regimes 



where the curvature is near the Planck scale (l p i = a/ Gh/c 3 w 1.6 x 10~ 33 cm). In 
a regime where the curvature approaches but always remains (significantly) less than the 
Planck scale, a semiclassical approximation to the full theory of quantum gravity should be 
sufficient. Examples of this regime include small evaporating black holes, when still much 
larger than the Planck mass (m p i = ^Tic/G 2.2 x 10~ 5 g), and the early universe after 
it has reached a size of many Planck lengths. In the standard semiclassical approximation, 
the gravitational field itself is treated classically, but is driven by the expectation value of 
quantum matter stress-energy. 

The form of the semiclassical corrections to Einstein's field equations is known for 
many important cases. 1 For example, for conformally flat classical backgrounds (in four 
dimensions), when the quantum state is constructed from the conformal vacuum, the 
corrections are completely determined by local geometry (the metric, the curvature, the 
covariant derivatives of the curvature): 2 

K(T ab ) = R ab - ^Rg a b + ^g ab 

+ a x h (^R 2 g ab - 2RRab - 2nRg ab + 2V a V b Rj 

(I 1 \ (L1) 

+ a 2 h ( -R cd R c d9ab - uR ab - -uRg ab + V a V b R - R cd R cad b j 

+ a 3 h (--^R 2 g ab + R cd R cadb ^j + 0(h 2 ). 

The parameters a\, «2, and as depend on the particular form of matter and regularization 
scheme, so we do not assume specific values or signs. Factors of % have been made explicit. 
Because the corrections are purely geometric, it is common to consider them not as matter 
source terms but as metric field terms (despite their matter origin). Non-conformally 
flat backgrounds can have more quantum corrections than Eq. (1.1). Of the examples 
above, the state-independent terms of Eq. (1.1) do not contribute in the case of the black 
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hole, where the exterior Ricci curvature vanishes, but they do contribute in the case of 
cosmological solutions. Several cosmological models are examined below. Because the new 
terms contain fourth derivatives in the metric, the new terms qualitatively change the field 
equations from system of second order equations to a system of fourth order equations. 

The new fourth order theory contains whole new classes of solutions unavailable to 
the classical theory. Many of these solutions have been examined. 3 One set of these so- 
lutions is particularly disturbing however. Solutions to the linearized theory around a 
flat background strongly indicate that flat space is unstable to ultraviolet fluctuations. 4 ' 5 
Using a 1/N approximation, Hartle and Horowitz showed that the ultraviolet instability 
can be made to occur at frequencies arbitrarily far below the Planck frequency, indicating 
that the instabilities cannot be easily fixed by calling the full quantum theory of gravity 
to the rescue. 5 Additional instabilities have also been found by Suen. 20 This strongly in- 
dicates that semiclassical gravity, if all its solutions are considered physical, is not a good 
description of the near classical limit of quantum gravity. 

It was shown in earlier work 6 that it is possible, and indeed desirable, to modify semi- 
classical gravity in a way suggested by and consistent with the perturbative nature of its 
derivation. The effective action and field equations of semiclassical gravity are perturba- 
tive expansions (formally, asymptotic expansions) in powers of h, truncated at first order 
in h. All behavior higher order and non-perturbative in fi has already been lost in the 
process of deriving the (approximate) effective action and field equations. Self-consistency 
then requires that only the solutions that are also asymptotic expansions in powers of K, 
truncated to first order, will be approximations to solutions of the full, non-perturbative 
effective action. Solutions not in this form are likely to be unphysical and should be ex- 
cluded. A simple model, presented below, will demonstrate that retaining non-perturbative 
solutions to a perturbatively derived higher derivative action results in false predictions. 
The non-perturbatively expandable solutions are spurious artifacts arising from the higher 
derivatives appearing in the perturbative correction, and will be referred to as spurious. 
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For convenience, perturbatively expandable solutions will sometimes be referred to as phys- 
ical, since only they correspond to predictions of the self-consistent semiclassical theory. 
For semiclassical gravity, it has been shown that the physical solutions show no signs of 
any instability of flat space (to first order in h). 6 

The easiest way of implementing the self-consistent method in semiclassical gravity is 
by reducing the fourth order equation, which has both physical and non-physical solutions, 
to a second order equation, which has only physical solutions (with one caveat described 
below). This iterative reduction has been demonstrated in a similar context by Bel and 
Sirousse-Zia for the case as = 0. 7 Much of the reduction (though not always all) can 
be done covariantly. It is clearly more efficient to find solutions to the reduced second 
order equations, almost all of which are physical, rather than finding all solutions to the 
full fourth order equations, most of which are spurious, and only using those which are 
physical. 

The aim of this work is to apply the reduction of order to a wide variety of gravitational 
systems. These include computing the reduced semiclassical equations for Friedmann cos- 
mologies (homogeneous isotropic solutions with perfect fluid matter), Friedmann-Lemaitre 
cosmologies (Friedmann cosmologies with cosmological constant), an interior Schwarzschild 
solution, and the general, conformally flat metric in terms of its conformal factor. Exam- 
ples of analytic and numerical methods are employed. In particular, we find the exact 
semiclassical solutions for spatially flat, radiation-filled Friedmann cosmologies, and exact 
and numerical semiclassical "bounce" solutions for radiation-filled Friedmann-Lemaitre 
cosmologies. As expected, the semiclassical corrections usually play only a small role in 
most systems far from the Planck scale. There are exceptions to this rule of thumb, how- 
ever, which we demonstrate by analyzing the semiclassical corrections to the (unstable) 
eternal Einstein universe. Here the corrections can cause large deviations from the classical 
solutions and yet remain within the domain of reliability. 
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We do not explicitly account for effects of particle creation (except in conformally flat 
spacetimes), only for "state-independent" contributions to the stress energy. 
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2. Review of Semiclassical Corrections 

The semiclassical field equations of general relativity (including cosmological constant) 
take the form 

Rab - ^Rg ab + Ag ab = K{T ah ) (2.1) 

where (T a b) = 0(h) is the expectation value or transition amplitude of the matter stress- 
energy tensor. For convenience, we consider only massless, conformally coupled fields (of 
arbitrary spin). We may reasonably restrict the form of (T &) to obey Wald's physical 
axioms: 8 

1) covariant conservation 

2) causality 

3) standard results for 'off-diagonal' matrix elements 

4) standard results in Minkowski space. 

Wald showed that any (T a &) that obeys the first 3 axioms is unique up to the addition 
of a local, conserved tensor. Furthermore, any local, conserved tensor can reasonably be 
considered part of the geometrical dynamics and so be written on the left-hand side of the 
field equations. We shall do so, rewriting Eq. (2.1) as: 

Rab - ^Rgab + Agab + Qab = K>{T ah ) (2.2) 

where Vt a b is conserved and purely local, i.e. it is constructed purely from the metric, the 
curvature, and a (finite number of) its covariant derivatives. 

Only terms in Q a b that are first order in % will be considered, since the semiclassical 
approximation already neglects higher order contributions. 21 Any term contributing to Q a b 
with a constant coefficient proportional to h must have dimensions of (length)~ 4 , since 
the only length scale is the Planck length, l p i, and % = l 2 x in units where G = 1. This 
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restricts the form of Vt ab for general spacetimes in four dimensions to linear combinations 
of two possible contributing terms. 9 



= \R 2 9ab - 2RR ab - 2nRg ab + 2V a V b R 



(2.3) 



{2) H ab = ^^- b jd^R^R cd 



(2.4) 



-R cd R cd g a b - uR ab - -Rg ab + V a V b R - R cd R cad b 



These two expressions are automatically conserved from their variational definitions. They 
are also fourth order in time derivatives of the metric. In conformally flat, four-dimensional 
spacetimes (where the Weyl tensor, C^e/> vanishes), ^H ab and ^H ab no longer remain 
linearly independent (in this case ^H ab = 3 ^ H ab ). However, a new quantity appears, 

^H ab = -—R 2 g ab + R cd R cadb 

12 (2.5) 

= —R^Rcb + \RR a b + \R C dR cd g a b — jR 2 g a b, 

which is conserved only in conformally flat spacetimes, but not as a result of a variational 
derivation, nor as the limit of a conserved quantity in non-conformally flat space-times. 10 
It is second order in derivatives of the metric, unlike ^H ab and ^H ab . Nevertheless, it 
is allowed by Wald's axioms, and, in general, contributes to the conformal anomaly. The 
most general expression for Vt ab is then 

Q ab = ai h {1) H ab + a 2 h (2) H ab + a 3 h {3) H ab + 0(h 2 ), (2.6) 

where it should be understood that the ^H ab term is only present when ^H ab = 3 ^H ab . 
Values of a±, «2, and as are predicted by specific matter couplings and regularization 
schemes, but we will treat them as free parameters. Factors of h have been made explicit. 
Inserting Eqs. (2.3)-(2.6) into Eq. (2.2) produces Eq. (1.1). 

As pointed out above, the new fourth order theory contains whole new classes of 
solutions unavailable to the classical, second order theory. It is the new solutions that 
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would indicate the instability of flat space. Since flat space is experimentally stable (or at 
least very metastable) , this strongly indicates that semiclassical gravity, if all its solutions 
are considered physical, is not a good description of the near classical limit of quantum 
gravity. However, the derivation of the theory is well founded, and it seems likely that some 
of the solutions do correspond to what we expect from quantum corrections to classical 
theory. It is necessary to break up the solutions to the semiclassical field equations into 
those we do not consider part of the theory (spurious or "unphysical" ) and the rest of the 
solutions ("physical"), which contain all the important information of the theory. 

That a theory should contain unphysical solutions should not be surprising to anyone 
who has examined Dirac's theory of charged particles including electromagnetic back- 
reaction. 11 The electromagnetic back-reaction problem shares several features with the 
quantum back-reaction described above. Its most important features are: 

1) A small correction term to the equations of motion changes the order of the equations 
of motion (from second order to third order) 

2) Not all solutions to the new equations of motion are physical; some must be excluded 
by external criteria. 

Dirac's equation of motion is 



The classic example of a non-physical solution has = but exponentially increasing 
acceleration, 



If the theory is to make useful predictions, unphysical solutions such as this must be 
excluded. In fact this so-called "runaway" solution shares much in common with the 
particular solutions to linearized semiclassical gravity which contribute to the instability 
of flat space. Other theories with higher derivative corrections, such as cosmic strings with 
rigidity corrections, can give unphysical solutions with negative kinetic energy. 12-14 
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(2.7) 




(2.8) 
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Similar problems occur even for semiclassical Quantum Electrodynamics (QED). The 
running of the electron charge coupling results in an effective action and Lagrangian density 
with higher order corrections 15 

I / Tie 2 \ 
£ eff = — A F^ [1 - ^^O) F lu , + ... + matter , (2.9) 

where e^ is the low energy electron charge. In this case it is clear what should be done. 
The higher derivatives in Eq. (2.9) do not correspond to new degrees of freedom for the 
electromagnetic field, but rather to the running of the charge. Since the effective action is 
a perturbative expansion in h, one may first solve the lower order equations of motion, and 
then solve higher orders iteratively. Treating Eq. (2.9) as a true fourth order equation would 
result in a theory very different from classical electrodynamics, where the electromagnetic 
field possesses negative energy modes. This is clearly undesirable from a physical point of 
view (though the resulting theory is mathematically well defined). 

One technique of distinguishing and excluding unphysical solutions of higher derivative 
theories, is called the self-consistent method, 13 ' 16 and was first applied in the case of the 
classical Dirac electron by Bhabha. 17 In the case of semiclassical QED shown above, it 
is equivalent to the obvious method of constructing higher order solutions iteratively. It 
can be applied to any theory for which higher derivative terms in the field equations are 
perturbative corrections to a lower order theory, and is most naturally applied to theories 
derived from an effective action with a small parameter that has been expanded in a power 
series. One expects extrema of that expanded action to be (perturbative) approximations 
to the extrema of the full effective action. For ordinary actions, whose variations give 
second order differential equations, this is usually true. For actions which possess higher 
derivative expansions, however, the opposite is true: most solutions of the perturbatively 
expanded field equations do not even have a perturbative expansion (i.e. they are not 
analytic in the expansion parameter as the parameter approaches zero). 
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An obvious cure for this behavior is to only rely on those solutions that are pertur- 
batively expandable (in the same sense as their effective action) to be physical. All other 
solutions are treated as spurious by-products of the higher derivatives, not to be considered 
part of the theory (e.g. the runaway Dirac electron). 

The self-consistent approach is extremely powerful. It removes all runaway and neg- 
ative energy solutions. It can be applied to semiclassical gravity as easily as to Dirac's 
classical electron. The amount of initial data required to specify a physical solution is the 
same as for the original uncorrected theory. 18 In the case of cosmic strings, where the full 
action is known exactly, any method not equivalent to the self-consistent approach simply 
gives the wrong results. The self-consistent approach seems clearly applicable to the case 
of semiclassical gravity. 

The solutions of semiclassical gravity are solutions to a fourth order differential equa- 
tion. The amount of initial data required to specify a physical solution, however, is the 
same as for classical gravity, which is given by solutions to a second order differential equa- 
tion. One can often find a second order differential equation which contains all the physical 
solutions to semiclassical gravity. This reduction of order greatly simplifies the process of 
finding physical solutions since most unphysical solutions are completely bypassed. The 
reduction is performed iteratively, using lowest (perturbative) order results to simplify the 
higher order semiclassical corrections. Reductions for corrections containing ^H a b and 
^H ab have been calculated for several cases by Bel and Sirousse-Zia. 7 The non-linearity 
of general relativity makes reduction of order awkward for the most general solutions un- 
less the stress-energy tensor has an extremely simple dependence on the metric. Just as 
in the case of classical gravity, however, the presence of symmetries can make soluble an 
otherwise intractable problem. We begin with the example of semiclassical corrections in 
Friedmann-Lemaitre universes. 

There is a small but important caveat regarding removing spurious solutions by the 
reduction of order. All the physical solutions to the higher order field equations are also 
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solutions to the reduced equations, but, still, not all solutions to the reduced equations 
may be physical solutions. If the corrections are nonlinear in the field variable, there may 
remain a much smaller number of spurious that must still be expunged. On reduction, 
however, it is often easy to identify the unphysical solutions and remove them. 

Reduction of order (and the above caveat) are well demonstrated by a simple example. 
The example uses a toy "semiclassical" equation 

x = -u 2 x + ah"x 2 + 0(ft 2 ), (2.10) 

which is non-linear in the fourth derivative correction. This is straightforwardly reduced 
with the substitution 

h'x = -hu 2 x + 0(h 2 ), (2.11) 

resulting in the reduced equation 

x = -uj 2 x + ahu 4 x 2 + 0(% 2 ). (2.12) 

This second order equation still contains all the physical solutions (perturbatively ex- 
pandable in h), but because of the quadratic nature of the equation, contains solutions 
non-perturbative in h. Solving the quadratic equation gives 

f - ~ 2 ^l = -oj 2 x + ahu 4 x 2 + OVh 2 ) 
x = 1+v ™^ . (2.13) 

The first of these manifestly produces only solutions perturbative in % and so contains 
only physical solutions. The second set is found only after treating Ti non-perturbatively: 
it would not be found using strictly perturbative methods and thus describes non-physical 
solutions. The point is that that even though Eq. (2.12) is reduced to second order, some 
care must still be exercised until the equation has been put into form of the first line of 
Eq. (2.13). In practice we will often use reduced equations in the form of Eq. (2.12), but 
if numerical methods are used, the form of Eq. (2.12) may not be adequate and one may 
need analogs of Eq. (2.13). 
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3. Friedmann-Lemaitre Models 
With Quantum Corrections 

In this section, we consider general Robertson-Walker metrics of the form 

ds 2 = -dt 2 + a(t) 2 (yY^2 + r2d ° 2 + r2 sin2 9d( t )2 ^j ' ( 3 - 1 ) 

where k takes the values ±1 or 0. The matter consists of radiation and the cosmological 
constant need not be zero. The (t, t) component of the Einstein equation with first order 
quantum corrections is of third-order in time derivatives of the scale factor a(t). We first 
reduce the equation to first-order (the same order as the corresponding classical Einstein 
equation), which the physically relevant solutions satisfy. For various values of the spatial 
curvature k and cosmological constant A, we obtain analytic and numerical solutions. 
For some solutions within this class of Friedmann-Lemaitre universes we find that there 
are models for which the first order quantum corrections remain small at all times. In 
some cases, the effect of small quantum corrections can cause a large deviation from the 
corresponding classical solution over a long period of time. 

In this conformally flat class of metrics, the Einstein equations with quantum correc- 
tions are of the form, 

R ab - (l/2)g ab R + Ag ab 

(3.2) 

+ ha x ^H ah + has {3) H ab + 0(h 2 ) = K T ab 

Here A is the cosmological constant, and T ab includes classical matter contributions and 
the lowest order expectation value of quantum matter, and ^H ab and ^H ab are defined 
bye Eqs. (2.3)and (2.5). The state-independent local quantum corrections of order % are 
included in the H terms on the left-hand-side. 
With the metric of Eq. (3.1), 

m „ -18A; 2 36/cd 2 54d 4 36d 2 d 18d 2 36da( 3 ) ^ ^ 

WHu = — + + — 3- + — — , (3-3) 

a 4 a 4 a 4 a 6 a z a z 
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and 



i^H tt = —r + —r- + —r. 3.4 

a 4 a 4 a 4 K 1 



The (t, £) component of the generalized Einstein equation (3.2), with matter consisting 
of radiation, is 

= -A-Kp ao/a 4 + (3/c + 3d 2 )/a 2 

+ aift(-18£; 2 + 36fcd 2 + 54d 4 - 36ad 2 d + 18a 2 d 2 - 36a 2 da (3) )/« 4 (3-5) 
+ a 3 h(3k 2 + 6ka 2 + 3d 4 )/a 4 + 0(h 2 ). 

Here, po is the radiation density at the time when the scale factor is ao- The a -4 dependence 
of the po term follows directly from conservation of the stress-energy and the fact that the 
the expectation value of the trace, T a a , vanishes to lowest order. The order h quantum 
corrections to this trace give rise to the ol\ and as terms. The values of of the ct's depend 
on which massless fields make up the radiation. In the spatially closed universe (k = 1), 
as a result of the non-local Casimer vacuum energy, the value of po will have added to it a 
small constant value proportional to h. 27 ' 28 There is no additional non-local contribution 
to the stress-energy because we are dealing with conformally invariant free radiation fields 
in Robertson- Walker metrics, which are all conformally flat. Before attempting to solve 
this equation, we reduce it to lower order in time derivatives of a(t), so that all solutions 
are assured to be of physical significance. 

3.1 Reduction of equation for scale factor 

Multiplying Eq. (3.5) by h and working to first order in h gives 

M 2 = -hk + h^- + h— + o(h 2 ). (3.6) 

3a 2 3 
This expression can be used twice to find that 

M* = hk 2 + h 2AKp : ao \ h K2po2a ° 8 



9 9a 4 
.2kKp a 4 ,2A/ca 2 A 2 a 4 2 
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(3.7) 



Differentiating Eq. (3.6) leads to 

ha = -h^^ + h^ + 0(h 2 ). (3.8) 

Here we have divided by a, assuming that it is not zero in the time interval of interest. 
Making use of this result twice gives 

ha = h ~ 9^ h 9a 2 + n_ 9~ + ° {fl } - (3 ' 9) 

We also need the third derivative, ha^ 3 \ which is obtained by differentiating Eq. (3.8): 

ha^ = h^ + + 0(h 2 ). (3.10) 

Using the above results repeatedly in Eq. (3.5) reduces it to the first order differential 
equation, 

3k 3d 2 KpociQ 4 8Ak/9 °o 4 



= -A+ — + — —. a x h 



a 2 a 2 a 4 a 4 ,„ ^ 

Except for the correction term involving a -8 , all of the terms of order % can be absorbed by 
a renormalization of the constants A and k. We make use of this renormalization when we 
discuss below models with radiation, spatial curvature and non-zero cosmological constant. 



14 



3.2 Spatially flat model with radiation and A = 

A simple but illuminating solution to the reduced Friedmann-Lemaitre-Einstein equa- 
tion (2.11) is the spatially flat (k = 0) case with zero cosmological constant (A = 0) and 
pure radiation (p = po(ao/a) 4 ). This should be a good approximation to our universe for 
one part of its history, after any inflationary epochs which smoothed out inhomogeneities 
but before massive fields cooled to non-relativistic temperatures. The classical solution, 
which has the scale factor grow as the square root of cosmological time, begins with a curva- 
ture singularity, and expands forever, becoming more and more flat. Because the physical 
semiclassical solutions are corrections to the classical solutions in powers of curvature, we 
expect that at late times, when the classical solution is nearly flat, the semiclassical cor- 
rections will be small. At early times, when the curvature is below the Planck scale but 
not above it, we expect the semiclassical corrections to be significant. At very early times, 
however, when the classical curvature is near or above the Planck scale, we expect that the 
semiclassical approximation will break down because neglected higher order corrections 
would dominate. We shall see how these effects manifest themselves below. 

For A = 0, k = 0, and p = po(ao/a) 4 : Eq. (2.11) reduces to 

which can be simplified by redefining the scale factor, cosmological time, and correction 
constant «3, in dimensionless units 



a = a 



i 



-l 



giving 



~ (AKp G al\ 2 _ / QAnpoa^ \ 



a = -^-ha 3 -^ + 0(h 2 ), (3.12) 



where a = This can be solved iteratively with the ansatz 



a(t) = a (t) + ha^t) + 0{h 2 ). (3.13) 
15 



Inserting Eq. (3.13) into Eq. (3.12) and expanding in powers of fi gives 

*o ~ 2 l~-2 
U : a = -a 

, , ' , (3-14) 
h : 2aoai = — -aia 3 + ct3a 6 . 

The first equation gives the one parameter family of (classical) solutions 

a (t) = (t-r )K 

where To is a constant of integration. When this is substituted into the second differential 
equation in (3.14), it gives the first order, linear, inhomogeneous equation 

5l + o7r~ — \ hl = "3(* - To)" S , 

2{t - T ) 

which, for a given tq, has a one parameter of family of solutions, 

_3 1 _i 
ai(t) = -a 3 (t - t ) 2 --ri(t-T ) 2 , 

given by T\, another constant of integration. The full solution is 

a(t) = (t-r )3 - ha 3 (t-T )-* - ^n(t - t )~^ + 0{% 2 ). (3.15) 

The apparently extra one parameter family of solutions (the freedom to specify T\ as 
well as To) is a result of the ambiguity in what it means to specify a solution. The 
freedom to specify both tq and T\ is equivalent to the freedom to specify one constant of 
integration, r, to lowest and first order in Ti. This can be seen by expanding r in powers 
of h, t = t + Titi + 0(h 2 ), and using it to write the solution (3.15) as 

a(t) = (t-r)* -ha 3 (i-T)-% +0(h 2 ). 



This freedom is related to the freedom of choosing different asymptotic expansions for the 
solutions, and is discussed in more detail in Appendix B. 
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Another solution to Eq. (3.12), to the same order in h, is 

A(t) = [(t - r) - 2ha 3 (t - r)- 1 ]? + 0(h 2 ). 

a and A differ only by terms higher order in h, so either solution is as good as the other. 
To the extent the two solutions disagree, the semiclassical approximation cannot predict 
which solution is more accurate. 

Fig. 1 shows plots of a, A, and do as functions of t, for two values of as (or, equivalently, 
as). The three regimes referred to above can be observed. When the universe is nearly flat 
(for i ~ 1), the semiclassical corrections are small. For intermediate scales (0.25 ~ t ~ 0.75, 
for the particular as plotted, as = ±0.01), the corrections are more substantial and cause 
noticeable deviations from the classical solution. For very early times (t ~ 0.25), the 
corrections dominate the classical solution. This, unfortunately, is the regime in which they 
cannot be trusted. This is made most obvious by noting the substantial differences (along 
the a scale) between the solutions a and A for small t. Since they are both equally valid 
solutions to the reduced semiclassical equations, to the extent they differ, the predictions 
of either are not meaningful. For positive (negative) as the effect of the corrections is to 
make the universe larger (smaller) than it would have been at small times, but the most 
dramatic predictions indicated by the solutions (and their plots) are in regimes where 
the solutions should not be trusted. All one can say is that at small times and large 
curvature, semiclassical corrections are important, and that at very small times and very 
large curvature, the contribution from higher order corrections of quantum gravity are 
necessary in order to make any meaningful predictions. The values chosen for as in Fig. 1 
(±0.01) have been chosen extraordinarily large (though still well within the perturbative 
regime) to demonstrate the qualitative effects of the quantum corrections. If one were to 
choose as of the Planck scale and a standard cosmo logical value for po^o; ^3 might be as 
small as 10" 100 . 
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3.3 Models with radiation, spatial curvature, and A^O 

More general models with general spatial curvature and cosmological constant and 
radiation, are governed to first order in % by Eq. (3.11). This can be rewritten in the form 

= ^ + ^- A ^-^T-+^ 3a8 +0(h), (3.16) 

where constant terms first order in h have been absorbed into a renormalization of the 
gravitational and cosmological constants: 

Ar = A(l - ^a 3 hA), (3.17) 



and 



2 

K r = k(1 a^h A r + 8aihA r ) . (3.18) 



With the exception of the final term, Eq. (3.16) has the same form as in the corresponding 
classical Einstein equation, but with k and A having constant corrections of order h. 

Before numerically integrating Eq. (3.16), it is advantageous to carry out some analytic 
simplification. Define the constants, 

A=^A ri (3.19) 

B = T^KrPoClo 4 , (3.20) 

C = BA~ l -Ak 2 A~ 2 . (3.21) 
Also define a dimensionless independent variable 

s = A 1/2 t, (3.22) 

and a function g(s) by 

g = a 2 -2kA-\ (3.23) 
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Then Eq. (3.16) can be written as 

&)W-c + »«, 4fa f 2 ff_ 1) , =a 0.24) 

Notice that C has the same value for k = ±1, so that the solution of this equation without 
the final term of order fi is the same for positive or negative spatial curvature, and is the 
same as the classical solution with renormalized gravitational and cosmological constants. 

Probably the most interesting case to consider in more detail is that of positive spatial 
curvature, k = 1, with positive cosmological constant, A > 0. This case includes both the 
"hesitation" and "turn-around" models of Lemaitre. 

The hesitation model spends a long time near the configuration of the Einstein static 
universe, during which density perturbations can grow rapidly. Close relatives of this 
model, which do not require a cosmological constant, have been recently considered in 
connection with the growth of perturbations and galaxy formation. 24 Such models may be 
affected by quantum corrections in a manner similar to that considered here. 

The turn-around model with no singularity is of interest because the perturbative 
approximation for the quantum correction is valid during the entire evolution of the model. 
Over a long period of time this quantum correction can cause a significant deviation of the 
radius of the universe from that of the corresponding classical model. 

The hesitation solution occurs when C > 0, and has zeroth order solution of Eq. 
(3.24), 

g = C 1/2 smh(s) (3.25) 
The turn-around solution occurs when C < 0, and has zeroth order solution, 

g = |C| 1/2 cosh(s) (3.26) 

In these solutions s could be replaced by s — s , with s constant. Also, since the zeroth 
order equation (for g) does not depend on the sign of k, these classical solutions remain 
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valid when k = — 1. However, the quantum correction does depend on the sign of k; and 
we will only consider the case of positive spatial curvature. 

We next consider the corrections to the turn-around model, and then the corrections 
to the hesitation model. 

3.3.1 Turn-around or bounce universe 

Dividing Eq. (3.24) by \C\, and defining the dimensionless function 

G(s) = |C|- 1/2 <7( S ), (3.27) 

and the dimensionless constants 

u = A-^Cl' 1 / 2 , (3.28) 
which regulates the abundance of radiation relative to the cosmological constant, and 

v = (l/4)ha 3 B 2 /(AC 2 ), (3.29) 

which regulates the quantum corrections, we obtain in the case of the turn-around solution 

(C<0): 

(^G(s)) 2 - G(s) 2 + 1 + v(G(s) + 2u)~ 2 = 0. (3.30) 
as 

(The hesitation model obeys the same equation with +1 replaced by —1.) The initial 
conditions for a turn-around solution are that G > 1 and -^G < 0, so that G will approach 
a minimum value near 1 and then increase. The perturbative solution will remain valid at 
all times for which 

v(G(s) + 2u)~ 2 < 1. (3.31) 
The relationship between a and G can be written as 

(A 1 / 2 a) 2 = (u- 1 G + 2) 1 / 2 . 
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The upper curves in Fig. 2 show the radius a as a function of time t, with both a and t 
measured in units of of A -1 / 2 = (4A r /3) -1 / 2 . Both curves have u = 0.5. The lower curve 
has v = 0. The upper curve includes the quantum corrections coming from 0:3, and has 
v = 0.4. The boundary condition for both curves was that G(l) = cosh(l). Again, as 
in Fig. 1, the numerical value of the quantum correction, in this case v = 0.4, has been 
chosen to be much larger than one would expect from Planck scale contributions, in order 
to emphasize the qualitative effects of the corrections. 

The curves were obtained by numerical integration of the equation for -^G(s). In 
solving Eq. (3.30) for -^G(s), the positive square-root was used in the range form s = 1 
to s = 0.09233 at which G(s) has its minimum value of 1.0467. The negative square-root 
was used for smaller values of s, since G(s) is a decreasing function of s in that range. The 
correction term in Eq. (3.30) remains smaller than 0.1 and has its largest value when G 
goes through its minimum. Thus, the perturbative solution is valid for all values of s. 

The lower curve is the same as the classical solution with A = A r , k = n r . It 
corresponds to the solution given in Eq. (3.26), and coincides with the quantum corrected 
solution at t = 1 (in units of A -1 / 2 ). It is clear that as a result of the deviation introduced 
near the minimum radius by the the quantum correction, the two models have increasingly 
different radii as one goes back into the past. 

3.3.2 Semiclassical Hesitation Universe 

From the previous examples, one can see that it is difficult to find examples of semiclas- 
sical corrections that are not overwhelmingly small. The last example used the exponential 
growth of the scale factor to magnify small correction terms. One instance semiclassical 
corrections might have a macroscopic impact is when the classical field equations describe 
an unstable system. Even initially small quantum effects could dominate the behavior at 
late times. 
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One important example of an unstable cosmological solution to the classical Einstein 
equations is the unstable Einstein static universe. Historically, Einstein looked for a static 
cosmological solution to describe an eternal universe, leading him to modify his equations 
with a cosmological constant. This solution is unstable, however, and small changes in 
initial data cause the universe to expand forever at an exponential rate (due to the cos- 
mological constant) or to collapse to a curvature singularity. One would expect quantum 
corrections to have similar, drastic effects. 

The simplest Einstein static universe is spatially closed (k = 1), contains a radiation 
fluid, with density p = po(ao/a) 4 , and cosmological constant, A, both carefully chosen 
such that A = 9/(4Kp ao)- Defining a natural length scale, L , such that L = a/3/A = 
a/4kpo«o/3, we can write the classical Einstein equation as 



0(h) = a 2 -±a 2 + k-^ 
3 3cr 

•2 t— 2 2 , i ' 

= a — L n a + 1 — 



(3.32) 



4Lo 2 a 2 



This has the solutions 



Ll/2 

a 2 = { {Ll/2) [1 - exp ( T 2(t - t )/L )] \ + 0{h), (3.33) 
{Ll/2) [l + exp(±2(f-f )/Lo)] 
plotted in Fig. 3. The first solution is the Einstein static universe. The second begins at 
a singularity and asymptotically approaches the Einstein static universe at late times (the 
positive exponential is the time- reversed solution) . The third spends an infinite time near 
the Einstein static universe, but pulls away and ends in an infinite inflationary epoch (the 
negative exponential is the time-reversed solution). The first and second solutions (and 
the time-reversed third) are unstable. 

The behavior of the semiclassical corrections should reflect the instability of the clas- 
sical solutions. The semiclassical Einstein equation of Eq. (2.11) becomes 

0(h) = a 2 - Ln 2 a 2 + 1 \t— - 6aift4r + a 3 hir ( L„ V H ^-^) . (3.34a) 

4L 2 a 2 a 2 a 2 \ 4L 2 a 2 / 
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This can be solved iteratively, as above, by expanding the solutions as a series in %: 

a = a c i + hai+0{Ti 2 ), (3.346) 



where a c i is one of the classical solutions in Eq. (3.33). Only the second and third solutions 
of Eq. (3.33) generate solutions to the semiclassical equations in general. [It should not be 
surprising that the constant solution has no semiclassical counterpart. The effect of the 
semiclassical corrections in Eq. (3.34a), if viewed as a Newtonian energy equation, is to 
change the potential term without changing the total energy. The energy corresponding 
to a static configuration without the corrections will not, in general, correspond to a static 
configuration when the potential is changed; to regain a static solution one must change the 
initial matter density appropriately] The equation of motion for at, obtained by inserting 
Eq. (3.34b) into Eq. (3.34a), is a linear, first order equation, and can be solved in closed 
form. 

The semiclassical solution generated by the initially singular classical solution is 
a 2 = ^(l-exp[-2(t-r)L- 1 ]) 



x 



l-exp[-2(t-r)Lo 1 ] 
. ^IhL- 2 (2 6Xp[2(t " r)L ° 1] exp[-2(t - r)Lo 1 



4 V l-exp[-2(t-r)Lo 1 ] (1 - exp[-2(t - t)Lq x ]) 2 (3.34) 

+ 8(t-r)V *M-2(t-T)L^ 
V ' l-exp[-2(t-r)Lo 1 ] 

+ 3 .""'f'r'u-'i ln(1 " exp[ " 2(t " 

1 - exp[-2(t - t)L ] /. 

+ 0(h 2 ), 

where r = to + ht\ + 0(Ti 2 ). As shown in Appendix A, the form of the solution (3.34) is 
ambiguous up to the addition of 0(h) terms proportional to d c z, arising from shifts 0(h) in 
the initial time (to — > r = t + hti). One may always choose a t\ such that the coefficient 
of the ambiguous term is zero. This solution is plotted in Fig. 4 for aiUL^ 2 = 0.0001 and 
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azhljQ 2 = —0.0001 (though the qualitative behavior is independent of those values, as long 
as 6«i — «3 > 0). As in the case of the spatially flat, radiation-filled universe above, the 
semiclassical approximation breaks down when too close to the initial singularity. At later 
times, where the semiclassical approximation is good, the effect of the corrections is to pull 
away from the Einstein static solution and begin an inflationary epoch. At very late times, 
when the correction terms dominate the classical solution completely, the corrections are 
untrustworthy, but at these late times the other (late-time-de Sitter) semiclassical solution 
exhibits the same inflationary behavior, and in a trustworthy regime. One can hope to 
match the two semiclassical solutions in the intermediate regime, where both are valid. 
We do this below. 

The semiclassical counterpart to the late-time-de Sitter classical solutions of Eq. (3.33) 
takes the form 



a 



2 L h + eM-^(t-r')L^]) 



x 



2 

2 exp^t-rO^o 1 ] 



1 — 3aihL 



l + exp[-2(t-r')Lo 1 ] 



+ —hLv 2 (2 6XP[2( * " T ' )L °" 1] exp[-2(t - t')Lo 1 



4 V l + exp[-2(t-r')L- 1 ] (1 + exp[-2(t - r')L^]) 2 (3.35) 

, ! exp[2(t - tQLq 1 ] 
( ] ° 1 + exp[2(t - r')L^] 

+ 3 i! xp[2 f 2 r /) ^ 1 - 11 ^ + ^-^ - ^o- 1 ]))" 

l + exp[2(t-r')iv ] /. 

+ o(h 2 ), 

where t' = t' -\-?it' 1 -\-O{?i 2 ) and t[ is chosen analogously to t\ above. Eq. (3.35) is plotted in 
Fig. 4 for two values of a.\ and ct 3 . At late times, the corrections are very small compared to 
the classical solution. At intermediate times, the corrections are small but non negligible, 
and at early times the corrections are so large as to be untrustworthy (ha\/a c i ~ 1). 

Because the semiclassical solutions of Eq. (3.34) and Eq. (3.35) are valid in different 
regimes, it is important to ask if there is any overlap of the regimes where both solutions are 
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valid. Furthermore, if there is such a regime, perhaps the solutions can be smoothly joined, 
corresponding to a universe beginning at large curvature near a singularity, flattening off 
at nearly constant scale factor for an extended period of time, and then proceeding to 
inflate in a de Sitter-like phase. This would correspond to a classical "hesitation" universe 
in which the matter density (or cosmological constant) is slightly greater than necessary 
for the Einstein static universe. 

For 6«i— a 3 > (the parameter range for which the solution of Eq. (3.34) is expanding 
at late times) there is an overlap region in which we can match the solution of Eq. (3.34) to 
the solution of Eq. (3.35), as shown in Fig. 4 by using the freedom to set the base times (r 
and t') of each solution individually. The matching can always be done smoothly, since the 
curves of a cross for all values of r — r', and we can adjust r — r' such that a is continuous 
(sufficient for matching solutions of a first order equation) . 

Furthermore, a is discontinuous only by terms 0(h 2 ). The matching can be done in 
regions where hai/a c i <C 1 for both solutions for a wide variety of parameters (such that 
6«i — as > 0). We may naturally interpret this joining of matched solutions as a unique 
solution to the semiclassical equation that is everywhere perturbatively valid (except the 
region near the initial singularity). The time of hesitation, th ~ t' — r, determined by 
the matching conditions, is logarithmically related to the coefficients of the semiclassical 
corrections: th ~ Lq In [(6cti — ot^KLiQ 4 ] , due to the exponential time dependence of the 
semiclassical solutions. 

The only potential obstacle to this interpretation is that, although Ua\ja c i <C 1 where 
the joining is done, ha\ja c i 1. We feel that this is no reason to doubt the validity of 
the joining, however, since ha\/a c i becomes large due to a c i vanishing, not due to hai 
becoming large. 
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4. Conformally Flat Corrections 



A conformally flat background spacetime (for which Eq.(3.2) gives the general form 
of the local corrections to the stress energy) has a metric tensor related to the Minkowski 
metric tensor by 



where e 2S is the conformal factor (a general scalar function), and i] ab is the Minkowski 
tensor (note that the conformal transformation is not, in general, a diffeomorphism) . Con- 
formally flat metrics are special in this sense — the entire metric (a symmetric tensor 
with 6 independent components) is completely specified by a single scalar function on the 
spacetime. 

As Eq. (4.1) is written, it is a tensor equation on the spacetime with classical metric. 
Here the Minkowski tensor r] ab is a tensor function defined on the physical spacetime. Since 
r\ a b is also the metric of flat Minkowski spacetime (or a piece of Minkowski spacetime), 
Eq. (4.1) implies that there is a map from the physical spacetime with metric g^ b to an 
unphysical flat spacetime with metric r\ ab , and the conformal factor e 2S can be viewed as a 
scalar function on either spacetime. If the conformal factor is known on the unphysical flat 
spacetime, this knowledge can be exploited to simplify the calculation of the semiclassical 
corrections, i.e. to do all the calculations on the unphysical flat background. 

We expand the general semiclassical Einstein equations for a conformally flat back- 
ground 



9ab = ^Vab, 



(4.1) 



G ab (g cd ) + a x h {1) H ab (g cd ) + a 3 h {3) H ab (g cd ) - KT ab (g cd ) = 0(h 2 ), 



(4.2) 



where 



9cd = 9 c cd + hh cd + 0(h 2 ). 



(4.3) 



This has the classical lowest order expansion 



G ab (g c a l b )-KT ab (g c J b ) = 0(h), 



(4.4) 
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which is already known if S is known. The equation first order in % is 

K S -^(9? d )hef + aih (1) H ab (gf d ) + a 3 h ^H ab (gf d ) - h^{gf d )h e f = 0(h 2 ) (4.5) 
og e f og ef 

where 

SG 11 1 

— —h ef = - -VaVbh - -Uh ab + V c V {a h b) c - -R cl h ab 

$9ef 2 2 2 (4 6) 

- \gt b [-Oh + V c V d h cd + R c J d h cd ] 

and all derivatives and raising of indices are with respect to the classical, physical metric. 

We have not explicitly expanded the last term of Eq. (4.5), the functional derivative of 

the stress-energy tensor, since its functional dependence on the metric depends on the 

particular form of matter present. If the functional dependence is known (as is often the 

case) then it is straightforward to calculate. 

This is a set of second order, linear, inhomogeneous equations for h ab . The second 
order equation produces a two-parameter family of solutions (just as for the classical equa- 
tion). The freedom to choose two additional free parameters in the semiclassical solutions 
arises from the freedom to specify the two parameters of the solution to the full semiclas- 
sical equations at both classical and semiclassical order independently. 

The power of Eq. (4.1) is most apparent when the conformal factor e 2S is known as 
a function of the fictitious flat spacetime. All semiclassical calculations can be performed 
on the flat spacetime instead of the physical spacetime. 

Then Eq. (4.6) becomes 

8G ab 1 _2s 



h e f = -^e [2/i c ( a; b) c — h- ab — h ab - c c + 8S- C c h ab + AS. c i a h b ^ 
<jg e f * 

— 4S' ; ( a /i b ) c; c + AS^ah-h) + 2S- C h ab . c 

+ 2S- a S- b h + S- c S- c h ab — 8S ]C S;( a h b ) c 

+ Vab(-h cd ., cd + h, c c + 2S, c h, c + 2S. c S, c h)l 



(4.7) 
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and the inhomogeneous terms are 

(1) H ab (9 c J d ) = {1) H ab (e^ Vcd ) 

= 6e~ 2S x [-12S. c S., a S ;b S^ - 12 S. a S ]b S ;c c + 12S. c S, b S. a c 

+ 12S, c S, a S, b c - AS, ca S, b c + 6S, b S, c c a + 6S, a S, c c b (4.8) 

— 4S- c S- a c b — 2S- C c ab + rj ab (3S- c S- d S^S d — 35 , ;c c 5 ;( j d 

- 12S, c s, d Sf d + AS, cd S^ d - 2S, c S^ d + 2S, C C /)] 

{3) H a M l d ) = ^H ab (e* s Vcd ) 

= e- 2S [-4S, c s, a S, b S, c - 4S, a S, b S, c c + AS, c c S, ab 

(4.9) 

+ 4S. c S. b S. a c + 4S. c S, a S. b c - 4S. ca S, b c 

+ Vab(S ;c S., d S, c S. d - 2S, c c S, d d - 4S, c S, d S, cd + 2S, cd S, cd )}. 
where semicolons refer to derivatives covariant with respect to the unphysical flat metric 
and all raising and lowering of indices in Eqs. (4.7)-(4.9) are with r\ ab . For a given stress 
energy, it is straightforward to put the last term of Eq. (4.6) into a similar form. By 
putting Eq. (4.6) in this form, the problem has now been simplified from solving a second 
order partial differential equation in curved spacetime to solving a second order partial 
differential equation in flat spacetime. 

Useful formulas analogous to Eqs. (4.8)-(4.9), but where the covariant derivatives and 
raising and lowering of indices are performed with the physical background metric have 
also been calculated. 19 These would be helpful in the case that the map from the physical 
spacetime to the unphysical spacetime implied by (4.1) were not known explicitly. 
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5. General Case 

The previous examples involved Friedmann-Robertson- Walker, Lemaitre, and other 
conformally flat spacetimes. We now turn to general spacetimes. The perturbative con- 
straints in the general case can be used to express the curvature tensors appearing in the 
first-order quantum corrections (i.e., the H : s) in terms of the lowest-order (classical) stress- 
energy tensor T ab of the matter. Usually this lowest-order T ab involves fewer derivatives 
of the metric than do the curvature tensors, so that this procedure results in an equation 
with fewer derivatives than the original. For example, for a classical fluid or a minimally 
coupled scalar field, the resulting equation contains no more than second derivatives of the 
metric. 

The reduction process must be modified when the stress-energy tensor contains explicit 
curvature terms, as for a conformally coupled scalar field. One way to deal with such a 
case is to evaluate the curvature tensor appearing in T ab using the lowest order classical 
solution, since only that will contribute to the correction terms of order h. 

For simplicity, we will suppose that a cosmological constant term, if present in the 
Einstein equation, is included in the definition of T ab . Such a term in T ab involves no 
derivatives of the metric. 

The perturbatively constrained equations should not have the instabilities exhibited 
when one tries to integrate numerically higher derivative equations. Such instabilities 
can be produced by the tendency for the growing or runaway solutions in the enlarged 
solution space to dominate the nearly classical perturbative part of solutions. The runaway 
solutions will not be present in the solution space after the reduction to lower derivatives. 

The first-order perturbative constraint coming from the Einstein equations with quan- 
tum corrections, Eq. (1.1), is 

hG ab = hKT ab + o(h 2 ), (5.1) 
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where T ab is the zero-order contribution of the stress-energy tensor. It follows that 

hR = -hKT + 0(h 2 ), (5.2) 



where T = T a a . Also, 

hR ab = hn(T ab - (l/2)g ab T) + 0(h 2 ). (5.3) 

Substituting these perturbative constraints into the expressions for the ^H ab gives (in 
four dimensions) the following results. 

^H ab = k (-2T. ab - (l/2)g abK T 2 + 2g ab T./ + 2 K TT ab ) + 0(h), (5.4) 
^H ab = k (-T. ab + T ab ./ + 2T a P. bp - (l/2) K g ab T 2 

(5.5) 

+ 2nTT ab - 2 K T pb T a p + (l/2) K g ab T pq T pq ) + 0(h), 

and 

^H ab = k (-(l/6)g abK T 2 + (l/3)nTT ab 

(5.6) 

- K T pb TJ + (l/2)K 9ab T pq T^) + 0(h). 
Then the perturbatively constrained Einstein equation with quantum corrections in a gen- 
eral spacetime is given by Eqs. (2.2)-(2.6) with these values for ^H ab and ^H ab (with 
A = 0). Only the lowest order or classical T ab appears in these quantum correction terms, 
and if there is a classical cosmological constant present, then it is included in the definition 
of the lowest order T ab . In general, the stress tensor contribution on the right- hand-side 
of the semiclassical Einstein equation will include non-local state-dependent contributions, 
such as those coming from gravitationally-induced particle creation and other effects. Local 
state-independent quantum corrections to the stress-energy tensor are, of course, already 
included with the local correction terms on the left-hand-side. If the zeroth order stress- 
tensor has vanishing trace T, as for radiation or massless particles, then the equations 
simplify considerably. 
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In a conformally flat spacetime, ^H ab is replaced by ( s >H ab , and further simplification 
may occur if the field is a massless conformally invariant field. The way in which this occurs 
was already discussed in the Introductory section. 

Similar expressions for WH ab and ^H ab were also obtained by Bel and Sirousse-Zia. 7 

5.1 Spherical Body With Quantum Corrections 

We next make use of the previous expressions for the correction terms, to write the 
equations governing the quantum corrections to the gravitational field of a static spherical 
body. These equations are in a form suitable for numerical integration.. 

We first calculate the local state-independent quantum corrections which enter into 
the Einstein equations for the most general spherically symmetric spacetime, which has 
the line element 

ds 2 = -B(r)dt 2 + A(r)dr 2 + r 2 d9 2 + r 2 sin 2 (6>)# 2 , (5.7) 

Here A and B will consist of a classical part and a quantum correction of order h. 

We will assume that the matter is described by a perfect fluid energy-momentum 
tensor, 

T ab = pg ab + ( p + p ) u a u \ (5.8) 

where u a is the four- velocity dx a /dr of the fluid volume element, and the proper pressure 
p and p are functions only of the radial coordinate r. Since the fluid is static, we have 
u r = u 9 = u* = 0, and u* = S(r)" 1 / 2 . 

From the previous section, we have the expressions for ^H ab: ^H ab , and ^'H ab 
in terms of the classical fluid T ab . Here we give the result of the calculation of those 
expressions (this is a lengthy calculation, but is simpler than calculating the quadratic 
curvature tensor expressions directly; in addition, because the perturbative constraints 
have been used, no higher than second derivatives appear in the result). 
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and 



One finds for the first correction term, 

^H rr = k [-(l/2) K p 2 A + npAp + (3/2)kA p 2 

+ 12p'/r + Sp'B'/B (5.9) 
-4p'/r -3p'B'/B] + 0(h), 

{1) H ee = k [-(l/2)Kr 2 p 2 + nr 2 pp + (3/2)«rV 
+ 6rp'/A - 3r 2 p'A'/A 2 

+ 3r 2 p'B'/(AB) + Qr 2 p"/A (5.10) 

- r 2 p'B'/(AB) - 2r 2 p"/A 

- 2rp'/A + r 2 p'A'/A 2 ] + 0(h), 

^H^ = sm 2 WHgo , (5.11) 

{1) H tt = [3KB/(2rA 2 )] [-nrp 2 A 2 + 2nrppA 2 
+ 3nrp 2 A 2 - 8p'A 

+ 2rp A' - Arp'A (5.12) 

+ Ap'A - rp'A' 

+ 2rp"A] + 0(h). 
The second correction term is found to be the following: 

^H rr = k [np 2 A + npAp - pB' 2 /B 2 

- pB' 2 /B 2 + 2p'/r + (3/2)p'B'/B (5.13) 

- (l/2)p'A'/A + p'B'/B + p"] + 0(h), 

( 2 > Hqq = k [nr 2 p 2 + kt 2 pp + rp' I A 

- (l/2)r 2 p'A'/A 2 + (l/2)r 2 p'B'/(AB) (5.14) 

+ rp '/ A + r 2 p"/A]+0(h), 
^H^ = sm 2 e ^H ee , (5.15) 
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and 

{2) H tt = k [3nBp 2 + ZKBpp - 2pB'/(rA) 

- 2pB'/(rA) + (l/2)pA'B'/A 2 

+ (l/2)pA'B'/A 2 - (l/2)pB' 2 /(AB) 

- (l/2)pB' 2 /(AB) + (l/2)p'B'/A (5.16) 
+ 2p'B/(rA) - (l/2)p'BA'/A 2 

- pB"/A - pB"/A 
+ p"B/A]+0(h). 

For the third correction term, which appears in place of ^ H ab in a conformally flat 
spacetime, the final result is much simpler: 

(3) # rr = (l/3)K 2 p(p + 2p)A + 0(h), (5.17) 

^H ee = (r 2 /A) ^H rr + 0(h), (5.18) 

WH 4>4> = sm 2 e i3) H ee , (5.19) 

and 

WH u = (l/3) K 2 p 2 B + 0(h), (5.20) 

It is interesting that ^'H ab is in the form of a perfect fluid energy-momentum tensor. 
It can be absorbed into the change, T ab — > T ab + AT ab , with 

AT ab = Apg ab + (Ap + Ap)u a u b + 0(h 2 ), (5.21) 

where u a is as before, 

Ap = -(l/3)a 3 hKp(p + 2p), (5.22) 

and 

Ap= -(l/3)a 3 hKp 2 . (5.23) 
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Therefore, in the case of a conformally fiat metric, the a 3 %^H ab quantum correction term 
in the Einstein equation, can be absorbed into a redefinition of the pressure and density: 

p -> p - (l/3)a 3 hnp(p + 2p) + 0(h 2 ), (5.24) 

and 

p -> p- (l/3)a 3 hKp 2 + 0(h 2 ). (5.25) 

The semiclassical Einstein equations with quantum corrections in general have the 
form, 

R ab - (l/2) 9ab R 

(5.26) 

+ «x h ^H ab + a 2 h ^H ab + a 3 h {3) H ab + 0(h 2 ) = K T ab , 
where T ab includes classical matter contributions and the lowest order state-dependent part 
of the expectation value of quantum matter fields. The state-independent local quantum 
corrections of order % are included in the H terms on the left-hand-side. It is understood 
that we may set a 3 = 0, except when the metric is conformally flat. In the latter case, 
it is understood that a<i = 0, since the first two corrections are then proportional to one 
another. In the conformally flat case, the a 3 term arises from the state-dependent part of 
the quantum stress-energy. 

With the expressions given above for the _£Ts, the Einstein equations are now easily 
written down for the general spherically symmetric metric. Only the state-dependent 
part of the expectation value of the quantum stress-energy tensor requires further work to 
calculate, but this will not increase the order of the highest metric derivative in most cases, 
so that the perturbative constraints have succeeded in reducing the semiclassical Einstein 
equations with quantum corrections to second order equations having the standard initial 
data. These equations are thus in suitable form for numerical integration. We will not 
carry that out here, but plan to return to it in a later paper. However, one spherically 
symmetric case where further simplification occurs will be discussed briefly in the next 
section. 
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5.2 Fluid Sphere of Constant Proper Classical Density 

Consider a fluid sphere which at the classical level has constant proper density. Let 
us suppose that, in addition to the classical fluid, only massless conformally invariant free 
fields, such as the photon and massless neutrino, are present. 

The classical interior solution for a fluid sphere of uniform proper density was found 
by Schwarzschild in 1916. 22 It is known that the Weyl tensor of this metric is zero, so that 
it is conformally flat. 

Because this spacetime is static and has no event horizons, we may suppose that 
the the quantum fields are in a well-defined vacuum state. It has been shown that for 
conformally invariant massless free fields in conformally flat spacetimes, the vacuum stress- 
energy tensor is determined by the trace anomaly. 23 ' 28 The vacuum stress-energy tensor of 
these fields is a linear combination of % ^H a i and h ^H a f,. We will suppose that there 
is no additional Casimir energy contribution in this spacetime. 

Therefore, the state-dependent part of the vacuum expectation value of the quantum 
stress tensor is zero for massless conformally invariant free fields propagating on this interior 
metric. The only effect of the quantum fields in their vacuum state is to give rise, through 
the conformal trace anomaly, to the ct\ and a 3 state-independent correction terms in the 
Einstein equations. Thus, in the interior of the fluid sphere, one has the complete equations 
which must be integrated. 

The classical interior Schwarzschild solution has the form of Eq. (5.7) with (for r < R) 

A cl (r) = (1 - 2GMr 2 /R 3 )~\ (5.27) 

and 

B cl (r) = ^ [3(1 - 2GM/R) 1 / 2 - (1 - 2GMr 2 /# 3 ) 1 / 2 ] . (5.28) 

The classical energy-momentum tensor corresponding to this solution is that of Eq. (5.8), 
with constant proper density, 

**> 
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The pressure p is 

(1 - 2GM/R) 1 / 2 - (1 - 2GMr 2 /R 3 ) 1 /' 



P{r) = 



(5.30) 



(1 - 2GMr 2 /R 3 y/ 2 - 3(1 - 2GM/R) 1 / 2 

For r > R, this interior metric joins with the classical Schwarzschild exterior solution of 
mass M, and zero density and pressure. It is known that of all stable fluid spheres having 
a given mass M and radius R, the Schwarzschild uniform density sphere has the smallest 
central pressure. For the pressure not to become infinite somewhere inside the object, it is 
necessary that GM < (4/9) R. This means that the radius of the static fluid sphere must 
be larger than the corresponding Schwarzschild black hole radius. Quantum corrections 
may possibly change the relationship between these two radii for sufficiently small fluid 
spheres. 

From the previous section, we have the expressions for ^H ab and ^H ab in terms of 
the classical fluid p and p. These are given for the constant density sphere by Eq. (5.29) 
and Eq. (5.30). The resulting semiclassical Einstein equations are of second order and are 
ready for numerical or analytic solution. We will carry this further in a later paper. 
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6. Conclusion 

We have considered first order semiclassical quantum corrections to a variety of clas- 
sical solutions to the Einstein gravitational field equations. We have used perturbative 
constraints to obtain the reduced semiclassical Einstein equations for Friedman-Robertson- 
Walker cosmologies, for Friedmann-Lemaitre cosmologies, for the gravitational fields of 
static spherically symmetric fluid bodies, and for the general, conformally flat metric in 
terms of its conformal factor. The reduced equations we obtained do not contain higher 
than second derivatives, and do not exhibit runaway solutions or instabilites of the original 
fourth order equations. They have the same physical content as the fourth order equa- 
tions, but yield only physically relevant solutions. Analytic and numerical solutions to 
these semiclassical equations were found in the cosmological cases. Although in most cases 
the semiclassical corrections play only a small role far from the Planck scale, there are 
some examples in which semiclassical quantum corrections cause significant deviation, or 
even qualitatively different behavior, from the classical solution. 

In the case of spatially flat, radiation-dominated Friedman- Robertson- Walker solu- 
tions, the corrections either strengthen or weaken the singular behavior at early times, in 
a regime where the perturbative corrections are valid (the perturbative validity does break 
down, however, before the time of the classical singularity itself can be reached). The cor- 
rections at late times become vanishingly small. In the Friedmann-Lemaitre "bounce" or 
"turn-around" solutions, quantum corrections can cause classical and semiclassical models 
which have the same initial conditions to have significantly different radii at late times. In 
the case of the "maximal hesitation" Einstein universe, the semiclassical corrections can 
cause large deviations and even qualtitatively different behavior from the corresponding 
classical solution. 

For fluid spheres, we have given the explicit first order quantum corrections as reduced 
field equations for the general case, and have discussed the constant classical density fluid 
in further detail. In future work, we intend to study solutions of these equations. For 
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small fluid spheres the corrections may significantly alter fundamental relations, such as 
the classical theorem which requires the radius of a static sphere of fluid to be larger than 
the radius of the Schwarzschild black hole having the same mass. 
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Appendix A 

Integration Factor 

In several instances in this paper, the semiclassical solutions are calculated iteratively, 
from an equation of motion of the form 

a(t) = f (a(t)) + %h{a{t)) + 0(ft 2 ), (Al) 

and an ansatz of 

a = a Q + hat + 0{% 2 ). (A2) 
Inserting Eq. (A2) into Eq. (Al) and expanding in powers of % produces 

d + hai = f {a ) + haxf^ao) + ft/i(a ) + 0(ft 2 ), 

which leads to the series of equations 

h° : ao(t) = f (ao(t)) 

h 1 : a 1 (t) = a 1 (t)ti(a (t)) + f 1 (a (t)). 
The first equation is typically a non-linear equation which might be solved in a variety of 
ways. The second is a first order linear inhomogeneous equation in a\{t) (once the classical 
solution ao has been determined), for which a general solution can be always found in the 
following form: 

ai(t) = ^T) + W)f dt ' MO/iMO), (A3) 

where c is an arbitrary constant of integration and (i(t) is an integrating factor given by 

/*(*) = «p(- f dt' /o>o(0))- (A4) 

Because Eq. (Al) has no explicit dependence on t, we know that there is a one parameter 
family of solutions, at (t) = a(t — to), to Eq. (Al), parametrized by the initial time to- The 
freedom to choose the constant c in Eq. (A3) must correspond to the freedom to change 
this initial time by to t = to + By making this shift, and expanding in powers 
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of h, we can determine the integrating factor without the need of integrating Eq. (A4) 
explicitly: 

ao(t - to) -> ao(t -t - Mi) = a (t - t ) - Mid (t - t ). (A5) 

Comparing this to Eq. (A3) reveals that [c//(t)] _1 = —tido(t — to), permitting us to rewrite 
Eq. (A3) as 

ai(t) = a (t) r^ /iMO) _ tlAo(t ) + Q ( fi ) (A6) 
J ao{V) 

where ti is an arbitrary parameter with dimensions of time, and the shift of initial time in 
Eq. (A5) induces only a higher order (0(h)) change in the integrand in Eq. (A6). 
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Appendix B 

Uniqueness of Perturbative Solutions 

For a general perturbative field equation of the form 

F (q, q, q) + eF 1 (q, q, q, q\ q iv ) + ... = 0(e n+1 ), (Bl) 

where e is the formal perturbative expansion parameter (e = h for semiclassical gravity) 
and q represents all the configuration space variables, there is some ambiguity in the way 
a perturbative solution, 

q = qo(t) +e qi (t) + ... + e n q n (t) + 0(e n+1 ), (B2) 

may be expanded in the same expansion parameter. For example, defining q as the 
quantity that satisfies the lowest order field equation, 

F (q ,q ,q ) = O(e) (B3) 

does not unambiguously determine qo, because we can shift by any quantity 0(e), i.e. qo — > 
qo + edq, and the new quantity will still satisfy Eq. (B3), only requiring an accompanying 
shift in the higher order terms in the expansion of the solution, i.e. qi — > qi — 5q. Similar 
ambiguity exists for the field equations of higher order terms 

fdF (q,q,q) dF (q, q, q) .. \ , iv . . 2 . 

v — dq — qi — 91 J + q,q > = °( e > 

etc. (B3a) 

Despite this ambiguity in breaking up the solution into terms that solve the field equations 
order by order, there is no ambiguity in the sum of all such terms. This can be seen by 
positing an additional requirement that the individual terms of the solution be explicitly 
independent of the perturbative expansion parameter, 

^ = 0( e " +1 -*). i = 0, ...,n (B4) 
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That this requirement can always be met can be easily seen as follows. Instead of solving 
Eq. (B3), solve the related equation 

F (qo,qo,qo) = 0(e n+1 ). (B5) 

Any solution to (B5) is also a solution to (B3), but there is no ambiguity to 0(e n ). Similarly 
solve the analogs of Eq. (B3a) to the highest order allowed 

fdF (q,q,q) dF (q,q,q). df (q,q,q)\ ...... lv n+1 

6 V dq qi M qi M 91 ) q,q ) = °( e ) 

etc. (B6) 

This process uniquely defines each of the terms in Eq. (B2), and therefore also the sum 
(to 0(e n )). 

Despite the ability to fix the expansion in this manner, it is often to our advantage 
to use the freedom to make order-by-order shifts in the terms qi of the solution, as done 
in Appendix A. Solutions with different e-dependencies can be obtained by adding e- 
dependent terms 

<?o — > go + tdiqo + e 2 5 2 qo + . . . + e n 5 n q 

qi^qi- $iqo - ... - e n ~ 1 5 n q + eSxqx + . . . + e n ~ 1 5 n - 1 q 1 

qi -> qi - hq\ - ... - e n ~ 1 5 n q 1 + €5^2 + • • • + e n ~ 1 5 n - 2 q2 

etc. (B7) 

where the 5\qj are arbitrary O(e ) functions to be chosen at one's convenience. 
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Appendix C 

Analytic Semiclassical Solution for Turn-around or Bounce Universe 

The semiclassical corrections to the turn-around or bounce universe were calculated 
numerically in section 3.3.1, but they can also be calculated analytically. The classical 
turn-around solution is a solution to Eq. (3.11), for positive cosmological constant (A > 0), 
spatially closed slicing (k = 1), and sufficiently small radiation density: 



< < 1. 

9 



(CI) 



The classical solution is 



1 3 

2 A 



+ 1 



where 



1 - 



9 



(C2) 



(C3) 



and < q < 1. 

The ansatz a = a c \ + Ua\ + 0(h 2 ), when inserted into Eq. (3.11) and expanded in 
powers of K, gives as the order % equation: 



0(h) = 2aia c i - 2—aia c i H — <na cl 



8aiAnpaQ _ 2 



a cl +a 3 



A 



2 , K P a -2 
~ a cl 



(C4) 



a 



This is a first order, linear, inhomogeneous equation in ai(t), and it may be solved by 
standard methods shown in Appendix A. The integrations are tedious, but easily within 
the grasp of a good symbolic integration software package. The result, up to the initial 
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time ambiguity dealt with in Appendix A, is: 



ai = Oil v6A((/ — q 1 ) cosh 
" 3 sinh 



2 4< 



gcosh 



2 4< 



-1/2 



+ 1 



+ «3- 



24 



2J$t 



Mn 



V 



a/1 — q 2 sinh 


2 \f 


3 1 


— cosh 


2^ 


- 9 


(/ cosh 


2 \ 


1 3 1 


+ 1 



A 
/ 



+ 2\/ -t sinh 



1 



H — 5ocosh 
9 



2./|t 



_1 __. 



Q q cosh 



x ! q cosh 

+o(n) 



2./ft 



-1/2 



+ 1 



+ 1 



(C5) 



The final result for a(t) is given by inserting Eqs. (C2) and (C5) into the perturbative 
ansatz, a = a c \ + fia\ + 0{Ti 2 ). 
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Figures 

Figure 1. Plots of the scale factor, a, A, and a as functions of t, for zero cosmological 
constant and zero spatial curvature. The top two plots are for as = +0.01, the bottom 
two are for as = —0.01, and the central plot, ciq , is for as = (classical solution), a and 
A are equally legitimate solutions to the semiclassical equations, differing only at order 
0(h 2 ). 

Figure 2. Plot of a as a function of t, each in units of (4A r /3) -1 / 2 . The lower curve is 
the classical solution (v = 0) with u = 0.5 and A = A r , k = n r . The upper curve includes 
quantum corrections, with v = 0.4 and other parameters unchanged. 

Figure 3. Solutions to the classical maximal hesitation equations. The upper plot is a 
maximal hesitation solution which spends an infinite time near the Einstein static universe, 
but pulls away and ends in an infinite inflationary epoch. The constant solution is the 
Einstein static universe. The lower plot, also a maximal hesitation solution, begins at a 
singularity and asymptotically approaches the Einstein static universe at late times. 

Figure 4. Solutions of quantum corrections to the maximal hesitation equations. Two clas- 
sical solutions and their semiclassical counterparts are shown for aiUL^ 2 = — asfiL^ 2 = 
0.0001. Both the classical and semiclassical initial singularity solutions begin at small 
scale factor, but the classical solution asymptotically approaches the static Einstein solu- 
tion while the semiclassical solution diverges exponentially from the classical solution (de- 
parting from the perturbative regime). The classical and semiclassical late time de Sitter 
solutions both agree at late times, but at early times the classical solution asymptotically 
nears the static Einstein solution, and the semiclassical solution diverges exponentially 
from the classical solution (departing from the perturbative regime). The two semiclassi- 
cal solutions may be smoothly matched, as seen in the insert, resulting in a non-maximal 
hesitation solution that is always perturbatively valid. 
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